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The effects of the variation of the dimensionless strong interaction parameter Xg — rUq/ Kqcd 
{niq is the quark mass, Aqcd is the QCD scale) are enhanced about 1.5 • 10^ times in the 7.6 eV 
"nuclear clock" transition between the ground and first excited states in the ^^^Th nucleus and 
about 1 • 10^ times in the relative shift of the 0.1 eV compound resonance in '^^"Sm. The best 
terrestrial limit on the temporal variation of the fundamental constants, \5Xq/Xq\ < 4 ■ 10~® at 
1.8 billion years ago {\Xq/Xq\ < 2.2 • 10"^^ y"^) , is obtained from the shift of this Sm resonance 
derived from the Oklo natural nuclear reactor data. The results for ^^^Th and '^^'^Sm are obtained by 
extrapolation from light nuclei where the many-body calculations can be performed more accurately. 
The errors produced by such extrapolation may be smaller than the errors of direct calculations in 
heavy nuclei. The extrapolation results are compared with the "direct" estimates obtained using the 
Walecka model. A number of numerical relations needed for the calculations of the variation effects 
in nuclear physics and atomic spectroscopy have been obtained: for the nuclear binding energy 
5E/E « —1.45 Smq/rriq, for the spin-orbit intervals 5Eso/Eso ~ —0.22 SrUq/mq, for the nuclear 
radius 5r /r « 0.3 5mq/mq (in units of Aqcd); for the shifts of nuclear resonances and weakly bound 
energy levels 5Er « 10 5Xq/Xq MeV. 
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INTRODUCTION 

Unification theories applied to cosmology suggest the 
possibility of variation of the fundamental constants in 
the expanding Universe (see, e.g., the review [ij). A re- 
view of recent results can be found, e.g., in Ref. In 
Ref. Q it was suggested that there may be a five orders 
of magnitude enhancement of the variation effects in the 
low-energy transition between the ground and the first 
excited states in the ^^^Th nucleus. This transition was 
suggested as a possible nuclear clock in Ref. Indeed, 
the transition is very narrow. The width of the excited 
state is estimated to be about lO""* Hz j5j]. The latest 
measurement of the transition energy gives 7.6 ± 0.5 
eV, compared to earlier values of 5.5 ±1 eV and 3.5 ±1 
eV Q. Therefore, this transition may be investigated 
using laser spectroscopy where the relative accuracy has 
already reached 10~^^. Several experimental groups have 
already started working on this possibility |^]. However, 
a recent paper [lo| claims that there is not any enhance- 
ment of the effects of the variation of the fundamental 
constants in this transition. The main aim of the present 
note is to demonstrate that the enhancement exists. We 
also estimate the relative shift of the 0.1 eV compound 
resonance in ^^"Sm to obtain new limits on the varia- 
tion of the fundamental constants from the Oklo natural 
nuclear reactor data [HI, [H, 13 1. 

We can measure only the variation of dimensionless pa- 
rameters which do not depend on which units we use. In 
the Standard Model, the two most important dimension- 
less parameters are the fine structure constant a = /fic 



and the ratio of the electroweak unification scale deter- 
mined by the Higgs vacuum expectation value (VEV) to 
the quantum chromodynamics (QCD) scale Aqcd (de- 
fined as the position of the Landau pole in the loga- 
rithm for the running strong coupling constant, as{r) ~ 
constant/ In (Agc_D?'/?ic)). The variation of the Higgs 
VEV leads to the variation of the fundamental masses 
which are proportional to the Higgs VEV. The present 
work considers mainly effects produced by the variation 
of Xq = rriq/ Aqcd where rriq = {jtIu + md)/2 is the aver- 
age light quark mass. Within Grand Unification Theories 
the relative variation of Xq may be 1-2 orders of mag- 
nitude larger than the variation of a [14]. Note that in 
the present work we do not consider effects of variation 
of the strange quark mass since they have larger uncer- 
tainty and should be treated separately. These effects 
were estimated in Refs. [sl [l5|. 



The results depend on the dimensionless parameter 
Xq — rriq/ Aqcd- In all calculations it is convenient to 
assume that Aqcd is constant and calculate the depen- 
dence on the small parameter niq. In other words, we 
measure all masses and energies in units of Aqcd and 
will simply restore Aqcd in the final results. Note that 
when a relative effect of the variation is enhanced it does 
not matter what units we use. The variation of the ratio 
of different units may be neglected anyway. 



2 



THORIUM 

To explain the origin of the enhancement we should 
present the small 7.6 eV interval between the ground 
and excited states in the ^^^Th nucleus as a sum of a 
few components which nearly cancel each other and have 
different dependence on the fundamental constants. If 
one performs the calculations exactly, it does not mat- 
ter how we select these components. However, in prac- 
tice the calculations are always approximate, therefore, 
a reasonable selection of the components will determine 
our final accuracy. For example, to study dependence 
on a we should separate the Coulomb energy from the 
remaining contributions to the energy. To study depen- 
dence on Xq = mq/AqcD it is convenient to separate out 
the spin-orbit interaction energy: 

u; = Eb + Eso^7.6eV . (1) 

Here Ei, is the difference in bulk binding energies of the 
excited and ground states (including kinetic and poten- 
tial energy but excluding the spin-orbit interaction) and 
Eso is the difference in the spin-orbit interaction ener- 
gies V/s(l • s) in the excited and ground states. We make 
this separation because we expect Eh and Eso to have 
a very different dependence on Xq = ruq/ Aqcd, as dis- 
cussed below. In ^^^Th the strength of the spin-orbit 
interaction is estimated to be Vis = —0.85 MeV from Ta- 
ble 5-1 of Ref. [16]. The difference of (1 ■ s) between the 
excited and ground states can be easily calculated using 
the expansion of the wave functions over Nilsson orbitals 
presented in Table 4 of Ref. 17]: Eso ~ 1.22V/s « -1-04 
MeV. (Note that without configuration mixing, for the 
"pure" Nilsson excited state [631]3/2+ and ground state 
[633]5/2+, Eso = 2V/S.) Then Eq. ^ gives us Eb « 
~Eso~l MeV and 



SoJ _ Esq ^ 6Eso _ 5Eb^^ - 1 3 iqS (^^12. _ 
CO LO Eso Eb Eso Eb 



(2) 



Qualitatively, we expect Eb and Eso to have a rather 
different dependence on Xj. In the Walecka model 
(which was used in Ref. [3| to estimate the enhance- 
ment factor) there is a significant cancellation between 
the a and u) meson contributions to the mean-field po- 
tential and the total binding energy E, while the a 
and w mesons contribute with equal sign to the spin- 
orbit interaction constant Vis fl8| . A similar argument 
may be made from the variational Monte Carlo (VMC) 
calculations with realistic interactions used in Ref. [19| 
to evaluate binding energy dependence on Xq. These 
calculations use nucleon-nucleon potentials that fit NN 
scattering data together with three-nucleon potentials 
that reproduce the binding energies of light nuclei. The 
binding energy is the result of a significant cancellation 
between intermediate-range attraction due to two-pion 
exchange and short-range repulsion arising from heavy 



vector-meson exchange. However, spin-orbit splitting be- 
tween nuclear levels has been found to be a coherent 
addition of short-range two-nucleon 1 • s interaction and 
multiplc-pion exchange between three or more nucleons 
[2Q] . Thus if meson masses move in the same direction 
due to an underlying quark mass shift, contributions from 
pion exchange and heavy vector-meson exchange will can- 
cel against each other in the binding energy, but reinforce 
each other in spin-orbit splittings. 



Binding energies 

The binding energy per nucleon and the spin-orbit in- 
teraction constant have a slow dependence on the nucleon 
number A. The total binding is dominated by the bulk 
terms, so we make the reasonable assumption that the 
variation of the bulk energy with Xq is the same for the 
two levels in ^^^Th and thus the variation of the difference 
SEb/Eb « SE/E. Moreover, the common factors (like 
jj^ tj^Q spin-orbit constant Vis [Si) cancel out in 
the relative variations SEso/Eso and SEb/Eb- Therefore, 
it may be plausible to extract these relative variations 
from th e ty pe of calculations in light nuclei performed 
in Ref. [l9|. The advantage of the light nuclei is that 
the calculations can be performed quite accurately, in- 
cluding different many-body effects. Their accuracy has 
been tested by comparison with the experimental data 
for the binding energies and by comparison of the results 
obtained using several sophisticated interactions (AV14, 
AV28, AVlS-hUIX - see (lil). As the first step, the varia- 
tions of the nuclear binding energies have been expressed 
in terms of the variations of nucleon. A, pion and vector- 
meson masses. The dependence of these masses on quark 
masses have been taken from Refs. 
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. The results 

for the relative variations of the total binding energies 
are presented in Tablc|T] (in the present work we add ^He, 
^He, and ^Be to this table). We see that all the results 
are close to the average value 5E/E « —1.45 SXq/Xq. 
The maximal deviations are for *He, which is especially 
tightly bound, and for ""He, which is a resonant state. 



Spin-orbit intervals 

To find the dependence of the spin-orbit constant Vis 
on rnq/h.QCD we calculate the spin-orbit splitting be- 
tween the pi/2 and levels in ^He, ^He, ^Li, and 
^Be in the present work. We use the Argonne vig two- 
nucleon and Urbana IX three-nucleon (AV18-I-UIX) in- 
teraction which provides our best results for small nuclei 
(see Ref. [3| for details and references). In all calcula- 
tions it is convenient to keep Aqcd = constant, i.e., mea- 
sure the quark mass niq in units of Aqcd- We restore 
^QCD in the final answers. As the first step we calcu- 
late the binding energies of the ground and excited states 
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TABLE L Dimensionless derivatives K 


_ SE/E 

SX^/Xq 


of the binding energy 


over Xq 


= mq/AgcD 








^He ''He '^He 


«He 


^Li ^He 


'Li 


'Be 


8Be 


^Be 


-1.39 


-1.44 -1.55 -1.08 -1.24 


-1.50 


-1.36 -1.93 


-1.50 


-1.57 


-1.35 


-1.59 



shown in Tabic HIl and their dependence on the nucleon, lowing energy shift for the 7.6 eV transition in ^^^Th: 
A, pion, and vector-meson masses, A£{niH) = — ^-^/-^ — 



8'mH I van ' 

shown in Table IIIII To find the dependence of these ener- 
gies on the quark mass, we utilize the results of a Dyson- 
Schwinger equation (DSE) study of sigma terms in light- 
quark hadrons |2l|. Equations (85-86) of that work give 
the rate of hadron mass variation as a function of the 
average light current-quark mass vriq — (ra^ -\- md)/1 as: 



5mH cfh 5mq 



rriH 



ruH ruq 



(3) 



with Giilnin values of 0.498 for the pion, 0.030 for the 
p-meson, 0.043 for the w-nieson, 0.064 for the nucleon, 
and 0.041 for the A. The values for the p and oj-mesons 
were reduced to 0.021 and 0.034 in subsequent work [23 |. 
We use an average of the p and lo terms of 0.030 for our 
short-range mass parameter my . 

It is convenient to present the result for the variation 
of the spin-orbit splitting in the following form: 



5E 



1/2 



5E- 



3/2 



E 



1/2' 



5Ei/-2 

'e^ 



E- 



5E- 



3/2 



3/2 



E- 



3/2 



(4) 



Accidentally, the calculated spin-orbit constant in ^He is 
the same as in 229Th, Vu = -0.83 MeV (the pi/2 - P3/2 
splitting in ^He is l.bVis). The spin-orbit constant in 
^Be is larger than in ^^^Th, in accord with the expected 
dependence A~^l'^ (see e.g. Ref. 
interval sensitivity coefficients Ks 



16|). The spin-orbit 



SEso 

E.^n 



i.so defined from 



rrin 



(5) 



for the quark mass variation in '^He, ^He, ''Li, and ^Be 
are -0.27, -0.16, -2.58, and -0.22, respectively. The 
^He, ^He, and ^Be values are all very similar, as all these 
nuclei are essentially one nucleon outside a 0+ core. The 
^Li value is anomalously large because its ground and 
first excited states are primarily a triton outside an al- 
pha core, so although dEso is comparable to ^Be, Ego 
is very small and not typical of the single-particle spin- 
orbit interaction we seek. Excluding the ^Li result gives 
us an average value of Kgo — —0.22 to use in ^^^Th. Note 
that the estimate based on the Walecka model, outlined 
in Sec. V below, gives a very similar value Kgo = —0.2. 



Frequency shift 

Substituting SEso/Eso = -0.22 SXg/Xq and 
SEb/Eb = -1.45 SXq/Xq into Eq. (O we obtain the fol- 



A V 

Soj = 1.2 i MeV 

Xn 



(6) 



This corresponds to the frequency shift Sv = 3 ■ 
10'^° SXq/Xq Hz. The width of this transition is 10""* Hz 
so one may hope to get the sensitivity to the variation 
of Xq about 10~^^ per year. This is 10^^ times better 
than the current atomic clock limit on the variation of 
Xq, '-^ lO"^** per year (see e.g. Ref. [3]). 



The corresponding relative energy shift is 
^ = 1.5.10^^. 

UJ Xa 



(7) 



This enhancement coefficient may be compared with the 
coefficient 0.4-105 from Refjl and 0.7-105 from Ref. [13]. 
The calculations in Ref. I23| have been done using the 
relativistic mean field theory (extended Walecka model) 
and some basic ideas from Ref. Thus, in this work we 
obtain an even larger enhancement! Here we present the 
relative variations from Refs. ^3i,i23.] for the new measured 
value 7.6 eV of the frequency u) (the old value was 3.5 
eV, and we multiplied the numbers from ^, '23| by (3.5 
eV)/w). The difference between the results of different 
approaches looks pretty large. However, this is only a 
refiection of the current accuracy of all three calculations. 
The present aim is to show that the enhancement does 
exist. 

Note that because of the huge enhancement it does not 
matter what units one will use to measure the frequency 
LO. In the calculations above we assumed that uj is mea- 
sured in units of Aqcd- However, the variation of the 
ratio of any popular frequency standard to Aqcd does 
not have such enhancement and may be neglected. 



Coulomb energy and effect of a variation 

We also would like to comment about the possible en- 
hancement of a variation. Ref. [13] claims that this en- 
hancement is impossible since the ground and excited 
states differ in the neutron state only and the neutron 
is neutral. Therefore, the ground and excited states 
have the same Coulomb energy and the interval does not 
change when a varies. We do not agree with this conclu- 
sion. Indeed, the total Coulomb energy of the ^^^Th nu- 
cleus is 900 MeV (see, e.g., [l^) which is 10^ time larger 
than the energy difference a'=7.6 eV. Therefore, to have 
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TABLE II: Experimental and calculated energies for the ground (P3/2) and first excited (pi/2) states of ^=5,7,9 nuclei in MeV. 

-•^Hed") ^He*(l-) 'He(|~) 'He*(^~) ^Li(a-) 'Li*(i~) "Bed') °Be*(i~) 

AV18+UIX -25.26 -24.02 -21.77 -19.56 -33.33 -33.02 -45.39 -42.01 

Expt. -27.41 -26.23 -28.83 -26.23 -39.24 -38.77 -58.16 -55.38 



TABLE III: Dimensionless derivatives A£(m_r/) = ^^^^^^ of the binding energy to the different hadron masses and the 
sensitivity K after folding in the DSE values of Smn /rnn ■ 







^Hed") 


'He*(i-) 


^Hed") 


'He*(i-) 




^Li*(i-) 


'Bed") 


''Be*(i~) 




+ 5n 


13.31 


13.83 


19.34 


21.34 


15.53 


15.48 


16.09 


17.12 


5a 




-10.24 


-10.72 


-14.92 


-16.63 


-11.96 


-11.88 


-12.39 


-13.27 




(+TNI) 


-5.82 


-6.07 


-8.78 


-9.73 


-6.91 


-6.88 


-7.27 


-7.76 


mv 




40.87 


42.84 


60.46 


67.54 


48.11 


47.81 


50.21 


53.85 


K = 


ii'/A' 

Srriq I m,q 


-1.24 


-1.29 


-1.93 


-2.13 


-1.50 


-1.49 


-1.59 


-1.70 



an enhancement it is enough to change the proton den- 
sity distribution (deformation parameter) by more than 
10~^. Any change in the neutron state influences the nu- 
clear mean field and proton distribution. For example, 
neutron removal changes the Coulomb energy of ^^^Th 
by 1.3 MeV This gives us an upper estimate (and 
a natural scale) for the change of the Coulomb energy in 
the 7.6 eV ^^^Th transition. One should expect a frac- 
tion of MeV change in any neutron transition in heavy 
nuclei. According to [l3| the weight of admixed octupole 
vibrations to the 7.6 eV state exceeds 20%. Octupole vi- 
brations involve both protons and neutrons. Therefore, 
the proton density distribution in the excited state is dif- 
ferent from the ground state and this difference is only an 
order of magnitude smaller than the difference in neutron 
distribution. 

The existence of the enhancement was confirmed by 
the direct calculation in Ref. [23] . The authors performed 
the calculation of the change of the nuclear mean field 
acting on neutrons induced by the change of the pro- 
ton density due to the variation of a. They obtained 
the enhancement coefficient 4 • 10"^. This corresponds to 
the Coulomb energy difference 0.03 MeV. As it was men- 
tioned above, an additional enhancement may come from 
the change of the nuclear deformation. Anyway, there is 
no doubt that the enhancement of the sensitivity to a 
variation in ^^^Th does exist. 



SHIFT OF THE RESONANCE IN SAMARIUM 
AND LIMITS FROM THE OKLO NATURAL 
NUCLEAR REACTOR DATA 



In Refs. 24 



|25| we derived a simple formula to esti- 
mate the shift of the resonance or weakly bound energy 
level due to the variation of the fundamental constants. 
Let us assume a Fermi gas model in a square well nuclear 
potential of the radius R and depth Vq. The energy of a 



single-particle energy level or resonance is determined as 

K,2 



2m 



(8) 



The momentum p in the square well is quantized, p 
constant/i?. Therefore, 



K 



2mi?2 



^0 



(9) 



For a resonance or a weakly bound level Er ^ { Er ■€i 
Vq) and the constant K « 2mR^Vo. Then we have 



SEr 



K Sm 2SR 
2mR^m^^' 



^^,5m 2SR SVo, 
m R Vq 
(10) 

This equation is also valid for a compound state with 
several excited particles. Indeed, the position of the com- 
pound state or resonance relative to the bottom of the 
potential well is determined mainly by the kinetic energy 
which scales as 1/R^ (both the Fermi energy and sum 
of the single-particle excitation energies scale this way). 
The shift of the resonance due to the residual interac- 
tion between excited particles (~ 0.1 MeV) is small in 
comparison with the depth of the potential well (Vq ~ 50 
MeV) and may be neglected. Note that the depth of the 
potential Vq is approximately the same in light and heavy 
nuclei. The radius of the well R w 1.2A^/^rQ, therefore, 
the relative variation SR/R = Sro/ro is the same too. 
Thus, the resulting shift of the resonance both in light 
and heavy nuclei is given by Eq. (jlOp and we may extrap- 
olate the accurate result for light nuclei to the resonance 
in isosm. 

In Table IIVI we present binding energies of the valence 
nucleon, S — —E (in MeV), and shift of the energy level 
(resonance), t— — = — x— — , due to the variation of 

the quark mass (in units MeV SXq/Xq) in light nuclei 
with A = 5,6, 7, 8, 9. In the derivation of Eq. pil)l it was 
assumed that the valence nucleon is localized inside the 
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TABLE IV: 

°He ^He ^Li '^He '^Li '^Be ^Be ''Be 

^expt -0.89 1.86 4.59 -0.43 7.25 5.61 18.90 1.67 

5'caic -2.24 -0.30 2.96 -3.19 5.11 3.52 16.89 -3.24 

- /•^°7* 3^42 9.92 9.52 11.73 15.35 15.53 17.24 16.19 

drriq /mq 

- /V" 1.62 6.12 7.06 4.58 11.62 11.45 15.82 6.52 



Based on the results above we conclude that l^rr^l < 

4-10^^ (for simplicity, we omit the small contribution of a 
variation here). Assuming linear time dependence during 
the last 1.8 billion years we obtain the best terrestrial 
limit on the variation of the fundamental constants 

1^1 <2.2.10-iV' • (16) 

A, 



potential well. This is not the case for ^He where the 
valence nucleon is localized mainly outside the narrow 
potential well produced by the *He core. As a result the 
potential < > averaged over the valence neutron wave 
function lps/2 is significantly smaller than the depth of 
the potential Vq. This explains why the shift in ^He (pro- 
portional to < y > - see Eq. (|10l) and Ref. [2^) is much 
smaller than the shift in other nuclei. Another extreme 
case is *Be where \Er\ is too large and the condition 
Er V is not fulfilled. The results for other nuclei are 
reasonably close to the average value 



SEr 



SX„ 
10 — 2- MeV 



(11) 



We assume this shift for the 0.1 eV resonance in ^^"Sm. 
This value does not contradict the order-of-magnitude 
estimates in Refs. 0, 24, [2^. Finally, we can add to 
this shift the contribution of a variation from Refs. f2fi| 
{5Er = -1.1 ± 0.1 MeV Sa/a). The total shift of the 
resonance in 



Sm is 



SEr 



10 (f^ 

Xq 



0.1 — ) MeV 

a 



(12) 



Now we can can extract limits on the variation of Xq 
from the measurements of SEr- Pioneering work in this 
area was done in Ref. We will use recent measure- 
ments [lll[l2,[l3| where the accuracy is higher. Ref. [l3| 
has given \SEr \ < 20 meV. Then Eq. ^ gives 



|^_0.1^|<2.10-^ 
X„ a 



(13) 



Ref. [IJ has given -73 < SEr < 62 meV. This gives 

|^_0.1^|<7.10-^ 
X„ a 



(14) 



Ref. [Ill has given -11.6 < SEr < 26.0 meV. This gives 

(15) 



, SX„ Sa , 
I— ^-0.1 — I < 2.6 -10" 
X„ a 



The limits on SEr have been been presented with 2a 
range. Note that Ref. jTlj has presented also the second, 
non-zero solution (it exists since the resonance has two 



tails): -101.9 < SEr < -79.6 meV. However, Ref. [13| 
tentatively ruled out this solution based on the data for 
the shift of a similar resonance in the Gd nucleus. 



VARIATION OF NUCLEAR RADIUS 

Variation of the nuclear radius is needed to calculate ef- 
fects of the fundamental constant variation in microwave 
atomic clocks where the transition frequency depends on 
a probability of the electron to be inside the nucleus. In- 
deed, the hyperfine interaction constant in heavy atoms 
has some sensitivity to the nuclear radius (including the 
Cs hyperfine transition which defines the unit of time, the 
second, and is used as a reference in numerous atomic and 
molecular clock experiments). This dependence was also 
requested by S. Schiller who proposed new experiments 
with hydrogen-like ions to search for the variation of the 
fundamental constants 27 1. 



In Table |V] we present a comparison of calculated and 
measured charge nuclear radii for the stable A = 2,7 
nuclei. Determination of the sensitivity of the nuclear 
radius to quark mass variation is a more involved calcu- 
lation than for the energy. While the dcuteron can be 
solved exactly, the VMC calculations for A > 3 nuclei 
of Ref. ^1^ have to be modified. This is because the 
variational bound for the energy is a quadratic function 
near its minimum in the space of variational parame- 
ters, but the radius is a linear function. In the previous 
VMC calculations, the variational parameters were fixed 
at the energy minimum for the nominal hadron masses 
corresponding to Smq = 0, and then not allowed to vary 
as the energy was evaluated for different Smn- Conse- 
quently the "size" of the trial wave function was essen- 
tially unchanged. For the radius determination, we must 
allow this size to vary. We do this by multiplying a set of 
variational parameters (those to which the radius is most 
sensitive) by a scale factor, and then carefully reminimize 
this scale factor for each Smn- This allows us to deter- 
mine Ar(TOff) — , — . The ASimn) — ^ ^^/^ — 
reported in Ref. |l9| are unchanged in this new mini- 
mization. 



TABLE V: Experimental and calculated point proton rms 
radii for stable A = 2 — 9 nuclei. 









=*He 


*He 


«He 


•^Li 


^Li 


^Be ^Be 


AV18+UIX 
Expt. 


1.967 
1.953 


1.58 
1.59 


1.77 
1.75 


1.45 
1.45 


1.92 
1.93 


2.46 
2.39 


2.34 
2.25 


2.45 2.40 
2.38 



This procedure works well for the A = 3, 4 nuclei, and 
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the Ar(mjy) is presented in Table IVll along with the to- 
tal sensitivity to the quark mass, Kr = ^^''/m ' obtained 
by folding in the DSE values for Smn/niH- However, be- 
cause our trial functions for A = 6 — 9 nuclei are not in- 
herently stable against breakup into subclusters, we need 
to make an additional constraint when calculating their 
sensitivity. For some of these nuclei, we have trial func- 
tions that asymptotically look like the appropriate sub- 
clusters bound in a Coulomb well with the experimental 
separation energy: ^Li is asymptotically an alpha and a 
deuteron bound by 1.47 MeV, ''Li is asymptotically an 
alpha and a triton bound by 2.47 MeV, and ^Be is asymp- 
totically an alpha and a ^He bound by 1.59 MeV. (^He 
is asymptotically a three-body a + n + n cluster and ^Be 
is an a -|- a -|- n cluster, so they cannot be treated this 
way.) For a quark mass shift Suiq/jriq ~ ±0.01, we know 
the total energy shift expected from our previous calcula- 
tions. We subtract that portion attributable to the alpha 
and deuteron or trinucleon subclusters, and use the re- 
maining energy shift to adjust the asymptotic separation 
energy of our trial function. This allows the size of both 
the subclusters and the well binding them to vary. For 
A = 6, 7, we have carried out this calculation for the total 
sensitivity only, and not for the individual Ar(TO//); 
these results are also given in Table IVIl 



TABLE VI: Dimensionless derivatives of point proton rms 
radii Ar(m£f) = ^^"^/^^ and the sensitivity with respect to 
mq after folding in the DSE values of Smu /niH- 









^He 


"He «Li 


■^Li ^Be 


rriN + 5n 


-7.32 


-4.81 


-4.73 


-3.04 




Sa 


4.07 


3.32 


3.28 


2.18 






2.57 


1.80 


1.77 


1.11 




my 


-16.39 


-12.97 


-12.79 


-8.50 




6r 1 r 
Sniq j rriq 


0.48 


0.34 


0.33 


0.20 0.35 


0.27 0.22 













The average value of Kr is about 0.3, which may serve 
as an estimate of the sensitivity for all nuclei. There are 
significant deviations from this value for the very weakly 
bound deuteron and very strongly bound "^He; the 
latter is probably a solid lower bound. 

The dependence of the nuclear radius on fundamen- 
tal constants manifests itself in microwave transitions in 
atomic clocks which are used to search for the variation 
of the fundamental constants (see e.g. Refs. 0, HSl)- The 
dependence of the hyperfine transition frequency on 
nuclear radius r in atoms with an external s-wave elec- 
tron is approximately given by the following expressions 
(in units of Kqcd)'- 



6r 



5m„ 



— Khr — KhrKr 

(jJh r nig 



0.3 K, 



hr 



nio 



Khr 



(27 - l)5h 
l-Sh 



(17) 



(18) 



(19) 



where 7 = (1 — Z^a^Y^^. For the Cs atom microwave 
standard the nuclear charge Z = 55 and Khr — —0.03; 
for the Hg+ microwave clock Z = 80 and Khr = —0.09. 

We also calculated the dependence of the "^He radius 
on a: j^j^ = 0.0034. For heavy nuclei the relative role 
of the Coulomb repulsion increases and the sensitivity to 
the a variation should be larger. 

ESTIMATES IN WALECKA MODEL 

It is instructive to compare the results obtained by 
the extrapolation from light nuclei with some "direct" 
calculations. In this section we estimate the variations of 
the resonance positions and spin-orbit splittings in heavy 
nuclei using the Walecka model [2^ where the strong 
nuclear potential is produced by scalar and vector meson 
exchanges: 



V 



2 -rms ^2 -rmv 



An 



An 



(20) 



Averaging Eq. (j20p over the nuclear volume we can find 



the depth of the potential well [24 1 
3 



Vo 



47rrj] 



9 s 



9v 



(21) 



Here 2ro = 2.4 fm is an internucleon distance. The result 
for the variation of the potential is 



Vo 



-7.5 



6ms 
ms 



5.5 



Smv „ Sro 



mv 



(22) 



Here we have used -%/-%- = 266.9/195.7 = 1.364 from 

Ref . [iBl ■ There is an order of magnitude enhancement of 
the meson mass variation contributions due to the can- 
cellation of the vector and scalar contributions in the de- 
nominator Vq . Eq. (|10p for the variation of the resonance 
position becomes 

<5£;^«y„(7.5^„5.5^-^ + ^). (23) 
ms my mN ro 

We do not know the variation of tq in the Walecka model, 
therefore, to make a rough numerical estimate we ne- 
glect this term. As above we take dependence of the nu- 
cleon and meson masses on the current light quark mass 
(m„ + ma) /2 from Refs [2l|, HI] : ' 



TO, 



-q 

nijv 



0.034 



Srrin 



0.064^, ^ = 0.013^, ^ = 0.498- 

JTlg m^ Trig ' m,y fftg 

The vector meson in the Walecka model is usually iden- 
tified with the w-meson so ^ = 0.034^. The scalar 
meson exchange, in fact, imitates both the a meson ex- 
change and two-pion exchange. Even if we neglect the 
two-pion exchange in zero approximation, there is vir- 
tual a decay to two tt. These virtual decays (loops on 
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a line in the NN-interaction diagrams with intermediate 
a) very strongly modify the a propagator and change its 
large distance asymptotics from e"'""'' to e"^™"'' [2^ . 
The mixing between a and two tt in mg should increase 
the sensitivity coefficient for the variation of mg. For an 
estimate we take an intermediate value between the neu- 
tron and vector meson mass sensitivity, ~ 0.05^^. 

ms m<j 

(Note that the positive contribution of ^ in Eq. 
produces an effect similar to that of an increase of ^^ifi-.) 
Then Eq. ^ gives 

- 10 — ^ MeV . (24) 

This rough estimate agrees with the result extrapolated 
from light nuclei. Note, however, that the accuracy of 
this estimate is very low due to the cancellations of dif- 
ferent terms. 

The scalar and vector mesons contribute with equal 
sign to the spin-orbit interaction constant Vis [18]. Also, 
the spin-orbit interaction is inversely proportional to the 
nucleon mass tojv squared. Thus, we have 

5Eso ^ ,SmN 0.58dms 0.42(5mv/ n o 
— = —2 ( I I ) « —0.2 . 

Ego rriN ms my ruq 

(26) 

This estimate is close to the result (—0.22 ^7^) obtained 
by the extrapolation from light nuclei. Note, however, 
that here we neglected the effect of variation of ro which 
probably should increase the absolute value of the sensi- 
tivity coefficient. 

CONCLUSION 

At the moment one can hardly calculate the sensitivity 
coefficient for the dependence of the strong interaction on 
the quark mass with an accuracy better than a fac- 
tor of 2. Moreover, it is hard to identify this dependence 
in phenomenological interactions which are used for the 
calculations in heavy nuclei. For example, it is not ob- 
vious that the scalar and vector mesons in the Walecka 
model are actually equivalent to free a and w mesons in 
particle physics. Therefore, to test conclusions obtained 
using the Walecka model, we explored a complementary 
approach. We performed the calculations in light nuclei 
where the interactions are well-known and the accuracy 
of the calculations is high. The binding energy per nu- 
cleon iJfc, the spin-orbit interaction constant Vis and the 
nuclear radius r have a slow dependence as a function of 
the nucleon number A. Moreover, the common factors 
(like A~^/^ in the spin-orbit constant Vis and A^/^ in the 
nuclear radius) cancel out in the relative variations 5r/r, 



5Vis/Vis and SEb/Eb. Therefore, we can extract these rel- 
ative variations from the calculations in light nuclei and 
use them in heavy nuclei. The errors produced by such 
extrapolation may be smaller than the errors of direct 
calculations in heavy nuclei. So far, this extrapolation 
and direct calculations using Walecka model give com- 
parable values of the enhancement factors in ^^^Th and 
isoSm. 
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